Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN M?eX style file vl.4) 



A ROSAT study of the cores of clusters of galaxies - I: 
Cooling flows in an X-ray flux-limited sample 

C. B. Peres, A.C. Fabian, A.C. Edge, S.W. Allen, R.M. Johnstone, and 

D. A. White 

Institute of Astronomy, Madingley Road, Cambridge CBS OHA 



00 

0\ 



1 February 2008 



> 

(N 
(N 

in 
o 

oo 

o 



ABSTRACT 

This is the first part of a study of the detailed X-ray properties of the cores of nearby 
clusters. We have used the flux-limited sample of 55 clusters of Edge et al. (1990) and 
archival and proprietary data from the iJO^S^T" observatory. In this paper an X-ray 
spatial analysis based on the surface-brightness-deprojection technique is applied to 
the clusters in the sample with the aim of studying their cooling flow properties. We 
determine the fraction of cooling flows in this sample to be 70-90 percent and estimate 
the contribution of the flow region to the cluster X-ray luminosity. We show that the 
luminosity within a strong cooling flow can account for up to 70 percent of a cluster X- 
ray bolometric luminosity. Our analysis indicates that about 40 percent of the clusters 
in the sample have flows depositing more than 100 yr~^ throughout the cooling 
region and that these possibly have been undisturbed for many Gyr, confirming that 
cooling flows are the natural state of cluster cores. New cooling flows in the sample 
are presented and previously ambiguous ones are clarified. We have constructed a 
catalogue of some intracluster medium properties for the clusters in this sample. The 
profiles of the mass deposited from cooling flows are analyzed and evidence presented 
for the existence of breaks in some of the profiles. Comparison is made to recent optical 
and radio data. We cross-correlate our sample with the Green Bank, NVSS and FIRST 
surveys and to the volume-limited sample of brightest cluster galaxies presented by 
Lauer and Postman (1994). Although weak trends exist, no strong correlation between 
optical magnitude or radio power of the brightest cluster galaxy and the strength of 
the fiow is found. 
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1 INTRODUCTION 

The X-rays emitted from clusters of galaxies represent a loss 
of energy of the intracluster gas. Thermal bremsstrahlung 
and line emission depend on the square of the gas density, 
which rises towards the cluster centre. As the density rises 
the loss rate increases and the resultant cooling timescale, 
tcooi decreases. Within the cooling radius (rcooi), the cooling 
time becomes smaller than the age of the cluster (tcooi < 
tagc ~ ^) and, under gravity and thermal pressure, the 
gas cools and flows inward to maintain pressure equilibrium. 

This simple picture describes the physics of a homoge- 
neous cooling flow. In a more realistic situation initial inho- 
mogeneities in the intracluster medium (ICM) lead to the 
establishment of a multiphase atmosphere which cools over 
a wide range of radii depositing cool material throughout 
the central ~200 kpc of a cluster of galaxies. (For reviews 
on cooling flows, see Fabian 1994, and Fabian, Nulsen & 
Canizares 1984). 



Although the mass cooling out of the X-ray band repre- 
sents only a negligible amount when compared to the total 
hot phase of the ICM, it is comparable to the mass of a 
brightest cluster galaxy (BCG) for flows which are undis- 
turbed for about a Hubble time. Therefore understanding 
the physical processes of a cooling flow and of the fate of 
the cooled material is of great importance for understand- 
ing the evolution of the core of a cluster. 

Edge et al. (1992) stress that the characterization of 
a cluster as a cooling flow depends heavily on the spatial 
resolution and sensitivity of the instruments used to image 
the cluster. This was the main limitation for early work. 

ROSAT offers, so far, the best opportunity to resolve 
the cores of the clusters and thus establish in a reliable man- 
ner the fraction of cooling flows in a flux-limited sample. It 
also offers a unique opportunity to study the occurrence 
of substructure in coohng flows (Peres, Buote & Fabian, 
in preparation) and to resolve mass deposition and cool- 
ing time profiles. The low background contamination of one 
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of its detectors, the Position Sensitive Proportional Counter 
(PSPC) allows surface brightness profiles to be extracted 
to a large radius. The good image capabilities of the other 
instrument, the High Resolution Imager (HRI), is impor- 
tant for the study of substructure and, in the context of the 
present work, fundamental for determining the properties of 
more distant flows. 

In the absence of a volume-limited sample, statistical 
studies of cooling flows are best pursued with the use of a 
fiux-limited sample, such as the one presented by Edge et 
al. (1990). The clusters in this sample are the brightest 55 
over the sky in the 2-10 keV band (B55 sample henceforth), 
and all have fluxes above 1.7 x 10~^^ ergcm^'^s"^. The 
clusters were selected from observations with the Einstein 
and EXOSAT observatories and the HEAO-1 and Ariel- V 
satellites. 

The work of Ebeling et al. (1996) on the X-ray Brightest 
Abell Cluster Sample (XBACS) showed that the complete- 
ness the B55 sample was satisfactory. A remaining issue in 
the analysis of cooling flow properties in this sample was the 
ability to image the cooling region in the centres of clusters. 
Here we analyse ROSAT pointed observations of the clus- 
ters in the B55 sample and deproject their surface brightness 
profiles. This allows us to build a catalogue of the ICM and 
cooling flow properties of the B55 clusters, and to compare 
these to recent optical and radio data. The structure of the 
paper is as follows. In sections 2 and 3 we describe the selec- 
tion of the observations and the method of analysis, respec- 
tively. In section 4 we present our results, and in section 5 
we summarize our conclusions. We assume throughout that 
Ho = 50 km s~^ Mpc^^ and go = 0.5. 



2 OBSERVATIONS 

We have used archival and proprietary observations of the 
B55 clusters carried out with the Position Sensitive Pro- 
portional Counter (PSPC) and the High Resolution Imager 
(HRI) on ROSAT. (The only exceptions were the cases of 
3C129 and A2147 for which no useful ROSAT observations 
exist ). The PSPC provides a FWHM spatial resolution of 
^ 25 arcsec (corresponding to a scale of 60 kpc for objects 
at a redshift z ~ 0.1) and a FWHM spectral resolution of 

AB 
E 



0.43 I 



, whereas the HRI has a FWHM 



0.93fceV; 

spatial resolution of ~ 4 arcsec and basically no spectral 
resolution. 

The HRI is ideal for characterizing mass deposition pro- 
files due to its high spatial resolution, whereas the PSPC is 
invaluable to recover ICM properties from the cluster sur- 
face brightness, due to its low and well characterized back- 
ground. In this work we have analysed 47 ROSAT PSPC 
observations, 38 ROSAT HRI observations (32 observations 
with both instruments), covering a total of 85 images of 53 
clusters. When more than one observation existed for the 
same object with the same instrument, we used only the 
observation with the longest exposure time. (Only in a few 
cases were other criteria, e.g. the existence of a better cen- 
tred image, considered more important.) 

A list of the clusters in the B55 sample, together with 
relevant generic information, is displayed in Tables ^ and ^ 
The X-ray positions given in Table Ol correspond to the peak 
in X-ray emission from our data, and were found by visual 



inspection. The optical positions were determined from op- 
tical counterparts in images retrieved from the Space Tele- 
scope Digitized Sky Survey (DSS). (We have used the HRI 
images, whenever available, to determine the locus of the 
X-ray peak.) For cases where no clear single X-ray peak ex- 
isted, like A1367 and A1736, the centroid of the X-ray emis- 
sion was used. The close agreement between the optical and 
X-ray positions ca n be seen from their offset, AS, presented 
in Table |l|(cf. sec. 

The radio information in Table |^ was obtained by 
searching the Green Bank 1400 and 4850MHz surveys (avail- 
able through NED) and the 1400MHz FIRST and NVSS on- 
line catalogues. The information on the existence of optical 
emission lines was obtained from a variety of sources in the 
literature and from private communication (C.S. Reynolds 
and C.S. Crawford). 

The X-ray observations used are identified by their 
Rosat Observation Request Sequence Number (ROR) in Ta- 
ble 1^, where the raw exposure times and cluster redshifts for 
each observation are also listed. The other entries in this ta- 
ble are discussed in the next section. 



3 METHOD OF ANALYSIS 

We have undertaken the present analysis in two distinct 
stages: (i) the reduction of the data to produce one image 
per observation, and (ii) the deprojection of the azimuthally 
averaged surface brightness profiles extracted from these im- 
ages. 

In the first stage, PSPC images have been constructed 
from counts in the energy range 0.4-2.0 keV and have been 
vignetting corrected. They have also been corrected for in- 
strument and telemetry dead time and had periods of high 
particle background and scattered solar X-rays removed. 
With the selection of this range of energies we minimize 
the Galactic and particle backgrounds (which dominate be- 
low 0.4 keV and above 2.1 keV, respectively), while re- 
taining most of the cluster emission. The vignetting and 
exposure-time corrections were performed using the IDL 
routines make.emap.pro, make jmage. pro, mexdiv.pro; this 
follows the prescription of Snowden et al. (1994). The back- 
ground subtraction of PSPC images was left to the deprojec- 
tion software, written by D. White and collaborators. (We 
select the bin for which the background counts are accu- 
mulated by inspection of the azimuthally averaged surface 
brightness, which works very well for the PSPC due to its 
low background count rate.) 

The HRI images were reduced in a different manner. 
Background subraction was performed with the STARLINK 
software ASTERIX, by accumulating counts in circles of ap- 
proximately 0.05 degrees in a clear off-centre area of the 
image. The vignetting correction was left to the deprojec- 
tion software. In all cases we were careful to inspect each 
background-subtracted surface-brightness profile and only 
deproject it up to a point where the difference between the 
subtracted and non-subtracted values differ insignificantly. 
In this way we guarantee that background subtraction does 
not affect the results significantly. 

For both the PSPC and HRI images we have identi- 
fied contaminating sources by eye and masked them out (we 
have also masked out the detector supporting structure in 
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Table 1. B55 sample, (a) Name of the cluster; (b) Right ascension of X-ray peak (J2000); (c) Declination of X-ray peak (J2000); (d) 
Right ascension of optical counterpart (J2000); (e) Declination of optical counterpart (J2000); (f) Offset between the X-ray and optical 
peaks, in arcsec; (g) Signal the presence {^), absence (x), or lack of information (— ) of optical emission lines; (h) Radio flux density 
(in mjy) at 6 cm from the Green Bank survey; (i) Radio flux density (in mjy) at 20 cm from the Green Bank (f), NVSS or FIRST (*) 
surveys. The values presented in columns (h) and (i) arc for sources which arc coincident with the optical or X-ray position within the 
angular resolution of the survey. We follow the convention used in (g) to indicate lack of information on radio fluxes. 



Cluster R.A. (X-ray) Dec. (X-ray) R.A. (opt.) Dec. (opt.) AO Opt. Lines Flux (5 GHz) Flux (1.4 GHz) 



(a) 






(b) 






(c) 




(d) 








(e) 


(f) 


(k) 


(h) 


(i) 


A85 


00 


41 


50, 


,8 


-09 


18 


07 


00 


41 


50, 


,4 


-09 


18 


12 


8 


V 


(4.6±l.l)xl0 


(5.8±0.3)xl0 


A119 


00 


56 


16. 


,8 


—01 


14 


45 


00 


56 


16, 


,1 


—01 


15 


19 


36 


X 




X 


A262 


01 


52 


45. 


,4 


36 


09 


26 


01 


52 


46, 


,5 


36 


09 


06 


24 


V 




(1.31)xl0^ t 


AWM7 


02 


54 


27. 


,4 


41 


34 


51 


02 


54 


27, 


,5 


41 


34 


46 


5 


X 






A399 


02 


57 


53. 


,6 


13 


01 


47 


02 


57 


53, 


,2 


13 


01 


50 


7 


X 




X 


A401 


02 


58 


56, 


,0 


13 


35 


03 


02 


58 


57, 


,8 


13 


34 


57 


27 


X 




X 


A3112 


03 


17 


57, 


,7 


-44 


14 


17 


03 


17 


57, 


,7 


-44 


14 


18 


1 


X 


(6.9±0.4)xl0^ 




A426 


03 


19 


48, 


,0 


41 


30 


46 


03 


19 


48, 


,3 


41 


30 


41 


6 


V 


(4.2±0.5)xl0* 


(2.12)xl0'' t 


2A 0335+096 


03 


38 


40, 


,2 


09 


58 


12 


03 


38 


40, 


,6 


09 


58 


11 


6 


V 




(2.41±0.12)xl0 


A3158 


03 


42 


50, 


,9 


-53 


37 


32 


03 


42 


53, 


,0 


-53 


37 


53 


28 


X 






A478 


04 


13 


25, 


,0 


10 


27 


59 


04 


13 


25, 


,3 


10 


27 


54 


7 


V 




(3.55±0.15)xl0 


A3266 


04 


31 


15. 


,7 


—61 


27 


08 


04 


31 


13, 


,5 


—61 


27 


12 


11 


X 






A496 


04 


33 


37. 


,6 


— 13 


15 


40 


04 


33 


37, 


,8 


— 13 


15 


43 


4 


V 


(4.4±l.l)xl0 




3C129 
















04 


48 


58, 


,2 


45 


02 


01 








(5.31) X 10^ 


A3391 


06 


26 


19, 


,9 


-53 


41 


53 


06 


26 


20, 


,4 


-53 


41 


36 


18 


X 


(1.922±0.019)xl0^ 




A576 


07 


21 


31, 


,2 


55 


45 


52 


07 


21 


30, 


,2 


55 


45 


40 


15 


X 




X 


PKS 0745-191 


07 


47 


30, 


,9 


-19 


17 


43 


07 


47 


31, 


,3 


-19 


17 


40 


6 


V 


(4.8±0.3)xl02 


(2.37±0.08)xl03 


A644 


08 


17 


25, 


,5 


-07 


30 


40 


08 


17 


25, 


,6 


-07 


30 


46 


6 


X 






A754 


09 


09 


18, 


,8 


-09 


41 


20 


09 


08 


32, 


,4 


-09 


37 


49 


691 


X 




(9.0±1.0)xl0" 


Hyd-A 


09 


18 


05, 


,8 


-12 


05 


40 


09 


18 


05, 


,6 


-12 


05 


44 


5 


V 


(1.40±0.01)xl0'* 


(4.08±0.13)xl0'* 


A1060 


10 


36 


43, 


,2 


-27 


31 


40 


10 


36 


42, 


,8 


-27 


31 


41 


6 


V 




X 


A1367 


11 


44 


48, 


,2 


19 


42 


05 


11 


44 


48, 


,0 


19 


41 


18 


47 


X 




X 


Virgo 


12 


30 


49, 


,0 


12 


23 


35 


12 


30 


49, 


,4 


12 


23 


26 


9 


V 


(6.1±0.8)xl0'' 


(2. 24) X 10* t 


Cent 


12 


48 


48, 


,9 


-41 


18 


44 


12 


48 


49, 


,1 


-41 


18 


42 


3 


V 


(1.53±0.08)xl03 




Coma 


12 


59 


35, 


,6 


27 


57 


31 


12 


59 


35, 


,6 


27 


57 


34 


6 


X 


(8.4±1.2)xl0 


(2.07±0.07)xl02 


A 1644 


12 


57 


12, 


,2 


-17 


24 


34 


12 


57 


11, 


,6 


-17 


24 


35 


9 


X 


(1.12±0.12)xl0^ 


(9.9±0.3)xl0 


A3532 


12 


57 


21, 


,8 


-30 


21 


51 


12 


57 


22, 


,0 


-30 


21 


50 


3 




(4.4±0.3)xl02 


(1.16±0.04)xl03 


A1650 


12 


58 


41, 


.7 


-01 


45 


44 


12 


58 


41, 


,5 


-01 


45 


41 


4 








A1651 


12 


59 


21, 


.7 


-04 


11 


47 


12 


59 


22, 


,5 


-04 


11 


46 


12 


X 






A1689 


13 


11 


29. 


,5 


-01 


20 


28 


13 


11 


29, 


,5 


-01 


20 


29 


1 


X 




X 


A1736 


13 


26 


50. 


,0 


-27 


10 


20 


13 


26 


48, 


,7 


-27 


08 


37 


103 


X 




X 


A3558 


13 


27 


56. 


,5 


-31 


29 


44 


13 


27 


56, 


,8 


-31 


29 


45 


4 


X 




(4.5±0.5)xl0° 


A3562 


13 


33 


36. 


,0 


—31 


40 


05 


13 


33 


34, 


,7 


—31 


40 


21 


23 


X 




X 


A3571 


13 


47 


28, 


,4 


-32 


51 


55 


13 


47 


28, 


,3 


-32 


51 


55 


1 






(8.4±1.7) 


A1795 


13 


48 


52, 


.7 


26 


35 


30 


13 


48 


52, 


,6 


26 


35 


35 


5 


V 


(2.6±0.3)xl02 


(9.3±0.3)xl02 * 


A2029 


15 


10 


55, 


,8 


05 


44 


46 


15 


10 


56, 


,1 


05 


44 


41 


7 


X 


(8.9±1.4)xl0 


(5.5) X 10^ t 


A2052 


15 


16 


43, 


,7 


07 


01 


19 


15 


16 


44, 


,6 


07 


01 


17 


13 


V 


(1.03±0.14)xl03 


(5.4) X 10-'' t 


MKW3s 


15 


21 


51, 


,8 


07 


42 


24 


15 


21 


51, 


,9 


07 


42 


30 


6 


V 




(1.26)xl0^ t 


A2065 


15 


22 


29. 


,0 


27 


42 


33 


15 


22 


29, 


,2 


27 


42 


26 


7 


X 





(1.4±0.3)xl0 


A2063 


15 


23 


04. 


,8 


08 


36 


20 


15 


23 


05, 


,3 


08 


36 


33 


15 


X 




(1.5±0.07)xl0 


A2142 


15 


58 


20. 


,2 


27 


13 


52 


15 


58 


20, 


,1 


27 


14 


00 


8 


V 




X 


A2147 
















16 


02 


17, 


,0 


15 


58 


27 




X 






A2199 


16 


28 


37, 


,7 


39 


33 


03 


16 


28 


38, 


,6 


39 


33 


04 


7 


V 


(4.8±0.5)xl02 


(3.7) xlO^ t 


A2204 


16 


32 


47, 


,1 


05 


34 


34 


16 


32 


46, 


,8 


05 


34 


31 


5 


V 




(7.01±0.2)xl0 


Tri Aust 


16 


38 


20, 


,3 


-64 


21 


28 


16 


38 


18, 


,3 


-64 


21 


36 


20 








A2244 


17 


02 


41, 


,9 


34 


03 


30 


17 


02 


42, 


,5 


34 


03 


35 


9 


X 




(2.41±0.13)xl0° * 


A2256 


17 


03 


13, 


,9 


78 


39 


06 


17 


04 


27, 


,1 


78 


38 


25 


59 


X 




X 


Ophiuchus 


17 


12 


27, 


,8 


-23 


22 


08 


17 


12 


28, 


,2 


-23 


22 


09 


6 


V 




(2.91±0.1)xl0 


A2255 


17 


12 


36, 


,2 


64 


04 


09 


17 


12 


35, 


,0 


64 


04 


14 


9 


X 




X 


A2319 


19 


21 


09. 


,7 


43 


56 


48 


19 


21 


10, 


,1 


43 


56 


43 


7 


X 




X 


Cyg-A 


19 


59 


28. 


,1 


40 


44 


05 


19 


59 


28, 


,4 


40 


44 


01 


5 


V 


(2.1±0.3)xl0S 




A3667 


20 


12 


24, 


,3 


-56 


49 


49 


20 


12 


27, 


,4 


-56 


49 


37 


25 


X 






A2597 


23 


25 


19, 


,3 


-12 


07 


20 


23 


25 


19, 


,7 


-12 


07 


27 


9 


V 


(4.1±0.2)xl02 


(1.88±0.06)xl03 


Klem44 


23 


47 


43, 


,4 


-28 


08 


20 


23 


47 


43, 


,4 


-28 


08 


37 


17 






(2.84±0.13)xl0 


A4059 


23 


57 


00, 


,2 


-34 


45 


39 


23 


57 


00, 


,5 


-34 


45 


35 


5 


V 


(1.17±0.13)xl02 


(1.29±0.04)xl03 
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Table 2. Input parametres for the deprojection analysis, (a) Name of the cluster. The letter inside brackets indicate wether 
the observation was made with the PSPC(P), or with the HRI{H); (b) Rosat Observation Request Sequence Number (ROR); (c) Raw 
exposure time in seconds; (d) Redshift of the cluster; (e) Bin size (in kpc) used in the deprojection; (f) Galactic neutral hydrogen column 
density along the line of sight to the cluster, given in units of 10^^ cm~^; (g) Temperature (in keV) from the catalogues of David et al. 
(1993) and White et al. (1997); (h) Outer radius to which the surface brightness profile was extracted (in Mpc). 



Cluster 


rxxJtx 




z 


bin size 




{kl ) 




(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


(h) 


A85(H) 


rh800271 


17308 


0.0521 


11.1 


0.3 


6.2 


0.244 


A85(P) 


rp800250 


10240 


0.0521 


41.5 


0.3 


6.2 


0.914 


A119(P) 


rp800251 


15203 


0.044 


35.6 


0.35 


5.1 


1.049 


A262(P) 


rp800254 


8719 


0.0164 


20.9 


0.53 


2.4 


0.452 


AWM7(H) 


wh800364 


14864 


0.0172 


3.8 


0.92 


3.6 


0.140 


AWM7(P) 


wp800168 


13335 


0.0172 


14.6 


0.92 


3.6 


0.532 


A399(H) 


rh800850n00 


6833 


0.0715 


73.5 


1.17 


5.8 


0.757 


A401(P) 


rp800235 


7465 


0.0748 


57.3 


1.11 


7.8 


1.175 


A3112(P) 


rp800302n00 


7600 


0.0746 


57.2 


0.4 


4.1 


0.886 


A3112(H) 


rh800627a01 


9540 


0.0746 


16.0 


0.4 


4.1 


0.421 


A426(H) 


wh800068 


10785 


0.0183 


5.93 


1.45 


5.5 


0.090 


A426(P) 


wp800186 


4787 


0.0183 


23.2 


1.45 


5.5 


0.766 


2A 0335+096 (H) 


rh800050 


14012 


0.0349 


5.14 


1.72 


3.0 


0.190 


2A 0335+096(P) 


wp800083 


10346 


0.0349 


28.7 


1.72 


3.0 


0.731 


A3158(P) 


rp800310 


3022 


0.0575 


68.1 


0.12 


5.5 


1.022 


A478(H) 


wh800091 


22712 


0.0882 


17.6 


1.36 


6.8 


0.520 


A478(P) 


wp800193 


22139 


0.0882 


66.0 


1.36 


6.8 


1.420 


A3266(P) 


wp800552n00 


13560 
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6.2 
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0.45 


5.2 
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rh800727n00 


10090 
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16.6 


0.56 


4.3 


0.236 
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wh800398n00 


23385 


0.1028 


20.0 


4.66 


8.5 


0.330 
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wp800623n00 


10477 


0.1028 


75.1 


4.66 


8.5 


1.464 
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rh800273n00 
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0.0704 


21.7 


0.73 


6.6 


0.355 
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wp800379n00 
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0.0704 


81.5 


0.73 


6.6 
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0.0124 
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rh800632n00 


14967 


0.0124 


11.0 


0.5 


3.3 


0.160 


A1367(P) 


rp800153 


18982 


0.0215 


45.2 


0.22 


3.5 


0.533 


VIRGO(H) 


wh700214 


14239 
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1.2 
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2.4 
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wp800187 


10539 


0.0037 
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0.25 


2.4 
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16526 


0.0109 
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0.8 


3.6 
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7985 


0.0109 


14.0 


0.8 


3.6 


0.331 
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rp800005 


22183 


0.0232 


29.2 


0.09 


8.0 


0.865 
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rh800242a04 


37410 


0.0232 


16.0 


0.09 


8.0 


0.250 


A1644(H) 


rh800851a01 


10232 


0.0474 


15.2 


0.47 


4.7 


0.137 


A3532(P) 


wp701155n00 


8620 


0.0585 


69.2 


0.62 


4.4 


0.876 


A1650(H) 


rh800852a01 


6054 


0.0845 


25.5 


0.15 


5.5 


0.297 


A1651(P) 


wp800353 


7435 


0.0846 


63.7 


0.17 


7.0 


1.243 


A1689(P) 


rp800248 


13957 


0.181 


115.6 


0.19 


10.1 


0.982 


A1689(H) 


rh800445n00 


13094 


0.181 


33.0 


0.19 


10.1 


0.823 


A1736(H) 


rh800853a01 


13783 


0.046 


29.6 


0.5 


4.6 


0.291 


A3558(H) 


wh800399 


16792 


0.0478 


25.6 


0.45 


6.5 


0.409 


A3558(P) 


wp800076 


30213 


0.0478 


96.0 


0.45 


6.5 


1.593 


A3562(P) 


rp800237n00 


20202 


0.0499 


59.9 


0.42 


3.8 


1.297 


A3571(H) 


rh800626n00 


19460 


0.0391 


25.5 


0.4 


7.6 


0.446 


A3571(P) 


rp800287 


6072 


0.0391 


31.9 


0.40 


7.6 


0.781 


A1795(P) 


rp80010-5n00 


36273 


0.0627 


49.1 


0.12 


5.1 


0.957 


A1795(H) 


rh800222a01 


11097 


0.0627 


13.0 


0.12 


5.1 


0.506 
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Table - continued. 



Cluster 


ROR 


At 


2 


bin size 




(kT> 


Rout 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


(h) 


A2029(H) 


rhl50024 


17757 


0.0767 


15.6 


0.31 


7.8 


0.383 


A2029(P) 


rp800249 


12550 


0.0767 


58.6 


0.31 


7.8 


1.025 


A2052{P) 


rp800275 


6215 


0.0348 


57.2 


0.29 


3.4 


0.557 


A2052(H) 


rh800223n00 


4429 


0.0348 


15.0 


0.29 


3.4 


0.308 


MKW3(P) 


rp800128 


9996 


0.0449 


72.5 


0.29 


3.0 


0.562 


MKW3(H) 


rh800425n00 


13502 


0.0449 


10.0 


0.29 


3.0 


0.244 


A2065{H) 


rh800724n00 


11780 


0.0722 


22.2 


0.29 


8.4 


0.363 


A2063{P) 


wp800184 


10198 


0.0350 


22.2 


0.29 


4.1 


0.363 


A2142(P) 


wpl50084 


7740 


0.0899 


67.1 


0.39 


11 


1.443 


A2142{H) 


rh800640n00 


19785 


0.0899 


27 


0.39 


11 


0.910 


A2199(H) 


wli800071 


5222 


0.030 


19.9 


0.09 


4.7 


0.195 


A2199(P) 


wp800644n00 


41082 


0.030 


24.8 


0.09 


4.7 


0.733 


A2204{H) 


rp800750n00 


15488 


0.1523 


40.9 


0.56 


9.0 


0.667 


A2204(P) 


rp800281 


5359 


0.1523 


102.2 


0.56 


9.0 


1.175 


TRI AUST(P) 


rp800280n00 


7338 


0.051 


40.7 


1.98 


7.9 


0.916 


A2244(P) 


rp800265n00 


2965 


0.1024 


112.2 


0.2 


7.1 


0.935 


A2256{P) 


wplOOllO 


17865 


0.0581 


45.8 


0.43 


7.5 


1.169 


A2256{H) 


rh800676n00 


43567 


0.0581 


25.0 


0.43 


7.5 


0.745 


OPHI(P) 


rp800279n00 


3932 


0.028 


23.3 


1.97 


9.0 


0.756 


OPHI(H) 


wli800067 


22285 


0.028 


9.3 


1.97 


9.0 


0.202 


A2255(P) 


rp800512n00 


14555 


0.0809 


92.0 


0.26 


7.3 


1.258 


A2319{H) 


wh800072 


5559 


0.0564 


29.7 


0.86 


9.9 


0.303 


A2319(P) 


wp800073a01 


3171 


0.0564 


44.6 


0.86 


9.9 


1.138 


C YG-A(Hj 


rn8UU021 


43240 


0.057 


24.0 


3.61 


7.3 


0.390 


CYG-A(P) 


wp800622n00 


9447 


0.057 


67.6 


3.61 


7.3 


1.464 


A3667{P) 


rp800234n00 


12560 


0.0530 


42.2 


0.4 


6.5 


1.328 


A2597{H) 


rh800111 


17996 


0.0824 


24.9 


0.25 


6.0 


0.241 


A2597{P) 


rp800112 


7243 


0.0824 


62.3 


0.25 


6.0 


0.904 


KLEM44(P) 


wp800354n00 


3353 


0.0283 


23.5 


0.15 


3.3 


0.505 


A4059{H) 


rh800224n00 


6320 


0.0478 


15.3 


0.11 


3.5 


0.190 


A4059{P) 


wp800175 


5514 


0.0478 


38.4 


0.11 


3.5 


0.595 



the PSPC images). Finally, we have extracted the profile of 
counts in annuli around the peak of the X-ray emission, to 
obtain a surface brightness profile suitable for deprojection. 
All the image analysis was carried out with ASTERIX and 
our deprojection software 

In the second stage we use the single-phase algorithm 
for surface brightness deprojection of Fabian et al. (1980), 
first used for the analysis of clusters of galaxies by Fabian 
et al. (1981)|^ In this technique counts in an X-ray image 
are eissumed to come from a spherically symmetric (single 
phase) hot thin plasma in hydrostatic equilibrium. Mod- 
elling the cluster emission as such, allows us to obtain the 
emission from each spherical shell in the cluster from the ob- 
served surface brightness profile and to compare this emis- 
sion with the predicted value from a plasma code. (For this 
work we have used the MEKAL code; cf. Mewe et al 1985, 
1986, Kaastra 1992, and Arnaud et al. 1985.) Through this 
comparison we can obtain the radial distribution of the ther- 
modynamic properties of the gas in the ICM, namely tem- 
perature, pressure and electron density. Once these quanti- 
ties are known, secondary properties like luminosity, central 
cooling time and mass deposition rates can be estimated. A 
more detailed exposition of the deprojection algorithm used 

* A more sophisticated deprojection, taking into account the 
multiphase nature of the ICM is known to give results in agree- 
ment with our simple approach here (Thomas et al. 1987). 



here can be found in White, Jones & Forman (1997) and 
references therein. 

The input parameters in the implementation of the de- 
projection algorithm used here are: redshift of the cluster, 
average temperature (as measured by Einstein MFC, EX- 
OS AT or GINGA; David et al. 1992), spectral band of the 
instrument (0.1-2.4 keV), Galactic hydrogen column density 
along the cluster line of sight, binsize of the data. Ho and qo, 
pressure at the outermost bin, cluster core radius and veloc- 
ity dispersion, galaxy linear mass (GLM), and galaxy linear 
mass cut-ofT (GCTOF) ^. From these, only four are kept free 
to fit the X-ray data: the cluster core radius and velocity dis- 
persion, GLM and GCTOF 0. We model the potential as an 

^ The gravitational model for the galaxy has a linearly increas- 
ing mass profile and a cut-ofi radius to avoid runaway. This is 
analogous to the model assumed by Thomas et al. (1987). 
■t Two caveats must be mentioned here: (1) Although we keep the 
value of the velocity dispersion free, the values used are roughly in 
agreement with (at least the lower limits of) the values measured 
by optical studies (cf. Zabludoff et al. 1990, Fadda et al. 1996). 
(2) The pressure at the outermost bin is not considered as a free 
parameter since it is always adjusted for the temperature at the 
outermost bin to agree roughly with the averaged value. (3) The 
parameters of the central cluster galaxy do not sensitively affect 
the results presented here in most cases, due to the resolution of 
the data. 
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isothermal sphere, adjusting the cluster core radius and ve- 
locity dispersion to conform to this model. This follows the 
works of Allen et al. (1996) and Fukazawa et al. (1994). The 
latter showed that the hot ICM phase exists up to the central 
parts of clusters of galaxies while the former showed that a 
consistent modeling of the cluster potential is achieved by 
assuming a roughly const ant average temperature. This is 
discussed at length in Sec. 



3.1 



Galactic absorption was accounted for by using the hy- 
drogen column densities quoted by Stark et al. (1992). No 
excess column density was introduced in our analysis, al- 
though this was found to exist from the X-ray spectroscopic 
studies of cooling flow clusters (White et al. 1991, Allen 
& Fabian 1997). Any excess absorption increases the mass 
deposition rates since the luminosity must be increased to 
produce the observed counts. 

We have usually rebinned the PSPC and HRI data to 
30 and 8 arcsec, respectively to minimize complexities asso- 
ciated with the instrument PSF (there are cases with coarser 
binning. The exact bin size used for each case can be found 
in Table Most of the values used for the average temper- 
ature were taken from David et al. (1992), although more 
recent data were used when available. The value used for 
A2597 is 6.0 keV. 

We have used a Monte-Carlo algorithm to perturb and 
deproject each surface brightness profile 100 times. By do- 
ing so we can place limits on the quantities derived from the 
deprojection analysis. However, we stress that these limits 
should be associated with the assumed potential only. Sys- 
tematic errors in the potential will lead to larger uncertain- 
ties in the derived quantities. 

As a final remark we note that the ICM is assumed to be 
in hydrostatic equilibrium. This is a key assumption needed 
to integrate the cooling fiow equations, once a gravitational 
potential is assumed. It is reasonable since the sound cross- 
ing time is usually shorter than the age of the cluster, which 
we assume to be ~ Hq^ ■ Supportive evidence for this and the 
hypothesis that the temperature profile is roughly constant 
is given below. 



3.1 The Multiphase (Inhomogeneous) Intracluster 
Medium 

Suggestive evidence that the highest temperature does not 
vary much with radius appeared in the work of Fukazawa 
et al. (1994), where the Centaurus cluster was observed 
with the ASCA satellite. They showed that hot components 
of the ICM exist throughout the central 5 arcmin (~ 130 
kpc) of this cluster. By fitting a two temperature Raymond- 
Smith(RS) model to the data these authors showed that the 
temperature profile for the cluster is approximately constant 
up to the inner ~ 2 arcmin (~50 kpc). A study of the cluster 
of galaxies A1060 by Tamura et al. (1996) reached similar 
conclusions. 

The hypothesis that the cluster temperature profile can 
be modeled as approximately flat received definitive support 
from the work of Allen et al. (1996). The authors studied the 
cluster PKS 0745-191 (0745 henceforth) with X-ray (spectral 
and imaging) and gravitational lensing methods. By carry- 
ing out a spatial analysis similar to ours the authors showed 
that the deprojected gravitational mass is in agreement with 
the mass inferred from the modelling of the lens provided 



the temperature profile is assumed to be fiat. In the study 
by Allen et al. (1996) additional supporting evidence for 
the assumption of a fiat temperature profile was obtained 
by spectroscopy. The X-ray spectrum of 0745, taken with 
ASCA, was only reasonably fitted if a two-temperature or 
a CFLOW model ^ were used. From this analysis Allen et 
al. (1996) showed that once the effects of the cooling flow 
are taken into account in the spectral analysis (through 
the CFLOW model) the cluster is observed to be approxi- 
mately isothermal. Comparing the X-ray and lensing results 
for 0745 the authors conclude that the temperature proflle 
should remain approximately isothermal up to ~10 arcsec 
(instead of the 2 arcmin limit obtained from the ASCA data 
alone) . 

From the work of Allen et al. (1996) it is clear that 
the assumption of an isothermal proflle for the gas makes 
the single-phase deprojection analysis appropriate for mod- 
elling the true multi-phase nature of the ICM and the cluster 
gravitational potential |^ (gravitational lensing was used as 
a cross-check) . For the clusters in this sample we cannot re- 
sort to the approaches discussed above to cross-check the 
appropriateness of our models, since spatially resolved spec- 
troscopic information is not available for most of the B55 
clusters, and their redshifts and masses are not suitable for 
gravitational lensing studies. However two indications of the 
appropriateness of our method are: (i) the size of the mass 
core radii we obtain for the B55 clusters and (ii) the func- 
tional form of the mass deposition rate profile inferred from 
the analysis. 

The values of mass core radii for the B55 clusters used 
in our analyses are available from Table 3 and are plotted in 
Fig- il- 

They were obtained from our requirement that the 
mean X-ray gas temperature stay approximatelly constant 
as the evidence presented above suggests. This methodol- 
ogy produces core radii which are overall smaller than ones 
found in previous analyses (White et al 1997, Edge et al. 
1992) where different assumptions were used. Our findings 
agree with the trend of smaller core radii found from grav- 
itational lensing studies of more distant clusters and with 
the multiphase picture of the ICM (Fabian et al. 1984; see 
also recent work of Miralda-Escude & Waxman 1995). 

Our model of the cluster gravitational potential in- 
cludes a separate mass component for the central galaxy. 
Despite this extra complication it is interesting to note from 
Fig. |l| that there is a clear separation between cooling-flow 
and non-cooling-flow clusters. This is supported by the X- 
ray/lensing analysis of Allen (1997). 



S This is the same model used in this paper and described in Sec. 
3. It is assumed that at each radial bin there is (i) gas in a single 
phase in hydrostatic equilibrium and (ii) gas cooling out of the 
flow. 

^ Allen et al. (1996) did not use a full multiphase deprojection 
method, resorting only to an additional cooling flow component 
(CFLOW) in their spectral analysis. This is exactly the compo- 
nent present in our single phase deprojection algorithm, and it 
simulates a whole range of phases cooling out of the X-ray band. 
Thomas et al. (1987) were able to show, through a study of 11 
clusters, that the single-phase deprojection method gives results 
in accordance with a multiphase treatment, where more than one 
phase (apart from the material which is dropping out of the flow) 
exists at each radial bin. 
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Table 3. Parameters for modeling the cluster potential. 

(a) cluster name, (b) cluster velocity dispersion (in kms~^), 
(c) cluster core radius (in kpc), (d) linear galaxy mass (in 10^ 
M0pc~^), (e) cut-of for the galaxy mass (in Mpc). Note that 
the profile for the Coma cluster was fitted with a single isother- 
mal model. 



200 400 600 BOC 

Core radius (kpc) 

Figure 1. The distribution of core radii in the B55 sample shows 
the same bi-modality found by the X-ray and lensing studies of 
Allen (1997). Cooling flows have smaller core radii overall and 
show a distribution peaked around 100 kpc. Non-coUing-flows dis- 
play a broader distribution of values, usually above 300 kpc. 



The widespread deposition of material over the cluster 
core also corroborates the multiphase nature of the ICM, 
where various different phases drop out of the flow at dif- 
ferent radii. In fact simple multiphase models proposed by 
Nulsen (1987) find M(< r) ~ r", with a ~ 1 (cf. Sec. 



4 RESULTS/DISCUSSION 

The deprojected profiles of all clusters analyzed by us can- 
not be presented here due to space limitations. Instead we 
present a list with the input parameters necessary to re- 
produce the analyses in Tables ^, and 0. The results are 
displayed in a catalogue form in Table M and discussed in 
the next subsections. 



4.1 X-ray Luminosity 

In Fig. ^ we plot the bolometric X-ray luminosity against 
total mass deposition rate. The bolometric luminosity was 
retrieved from the work of David et al. (1992) and the mass 
deposition rates come from our analysis. We have used HRI 
observations here only when no PSPC observation was avail- 
able. From this plot we notice that the upper left and the 
lower right corners are not populated. The former can be 
taken as an expression of the fact that for a given X-ray 
luminosity there is always a maximal mass deposition rate. 
(Note that the cooling luminosity Lcooi ~ M x T.) The lat- 
ter gives us a clear indication that very luminous systems 
harbour massive flows. 

We have also investigated how much of the total bolo- 
metric luminosity can be attributed to the flow for each 
cluster in the sample. To answer this question we computed 

the luminosity from within the cooling radius, L{< rcooi) f\ 
and plotted it against the bolometric X-ray luminosity given 



Cluster 


'-''dcproj 


'^dcproj 


GLM 


GOTO 


A85 


700 


130 


2.0 


0.06 


A119 


680 


600 


1.0 


0.08 


A262 


450 


130 


1.2 


0.04 


AWM7 


520 


200 


1.1 


0.08 


A399 


940 


550 


1.0 


0.10 


A401 


950 


350 


0.8 


0.08 


A3112 


670 


120 


1.4 


0.06 


A426 


700 


60 


0.9 


0.06 


2A 0335-1-096 


600 


40 


0.1 


0.05 


A3158 


860 


400 


0.8 


0.08 


A478 


850 


180 


1.8 


0.08 


A3266 


780 


500 


1.5 


0.17 


A496 


610 


60 


1.2 


0.05 


A3391 


750 


300 


0.8 


0.08 


A576 


700 


150 


0.6 


0.06 


PKS 0745-191 


1050 


100 


1.5 


0.08 


A644 


860 


200 


0.6 


0.08 


A754 


800 


300 


1.6 


0.09 


HYD-A 


680 


220 


1.8 


0.07 


A1060 


500 


150 


0.7 


0.09 


A1367 


700 


500 


0.4 


0.03 


VIRGO 


600 


90 


1.0 


0.02 


CENT 


400 


120 


1.6 


0.06 


COMA 


1000 


500 


- 


- 


A1644 


700 


180 


1.8 


0.08 


A3532 


660 


350 


0.8 


0.12 


A1650 


600 


100 


0.8 


0.08 


A1651 


900 


220 


1.0 


0.12 


A1689 


1150 


200 


1.4 


0.15 


A1736 


800 


700 


1.2 


0.05 


A3558 


970 


600 


2.0 


0.15 


A3562 


550 


180 


0.8 


0.12 


A3571 


850 


220 


1.1 


0.08 


A1795 


710 


70 


0.5 


0.07 


A2029 


920 


100 


0.8 


0.06 


A2052 


550 


200 


2.1 


0.08 


MKW3s 


550 


100 


0.8 


0.08 


A2065 


750 


200 


2.1 


0.10 


A2063 


600 


150 


1.0 


0.08 


A2142 


1050 


200 


1.8 


0.12 


A2199 


670 


150 


1.3 


0.08 


A2204 


1100 


80 


1.0 


0.08 


TRI AUST 


800 


220 


2.5 


0.05 


A2244 


930 


150 


0.5 


0.05 


A2256 


920 


550 


0.9 


0.15 


OPHI 


880 


220 


2.5 


0.06 


A2255 


1000 


650 


0.5 


0.07 


A2319 


870 


180 


0.4 


0.08 


CYG-A 


850 


50 


2.5 


0.05 


A3667 


730 


250 


0.5 


0.04 


A2597 


850 


100 


1.6 


0.10 


KLEM44 


550 


100 


0.6 


0.08 


A4059 


560 


110 


0.4 


0.09 



Note that L{< r^^ol) includes direct cooling and the gravi- 
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Figure 2. Mass deposition rate vs. bolometric X-ray luminosity for the clusters in our sample. The lack of points in the upper left 
corner only reflects the maximum M allowed for a given luminosity. Stars are used to denote PSPC observations. Squares represent the 
clusters for which only HRI observations were available. No special symbols are used to denote the clusters which have mass deposition 
rates consistent with zero. The error bars for these clusters extend to the bottom of the figure and are kept to indicate the upper limit 
in their mass deposition rates (When the median mass deposition is zero we assume it equal to 0.1 for display purposes). The luminosity 
of the Virgo cluster was taken from White ct al. (1997) . 



by David et al. (1992) in Fig. |. We find that, despite tiie 
spread, 20 percent of tiie clusters have 50 percent or more 
of their bolometric X-ray luminosity produced in a region 
< ^cooi- The cooling flow region in A2204 contributes more 
than 70 percent to the cluster's total X-ray luminosity. The 
percentages for the other clusters are listed in Table o. 

Tfie approach described in the last paragraph does not 
treat non-cooling-fiows properly since the concept of a cool- 
ing radius, rcooi, in this case is meaningless. To compare 
cooling flows and non-cooling-flow clusters we have com- 
puted the fraction of the total luminosity coming from a 
central region of 100 kpc radius. This quantity, /loo, is plot- 
ted against M in Fig. ^ from which is clear that: (i) non- 
cooling-flows have an overall smaller fraction of their X-ray 
luminosity emitted from the centre (/loo ~ 1 — 10 percent), 
when compared to cooling-flows (/loo ~ 5 — 50 percent; clus- 
ters with massive cooling flows show /loo > 15 percent) and 
(ii) a trend of increasing /loo with M is indicated by the 
plot Some of these issues were investigated previously by 



tational work done on the gas; direct cooling accounts for only 
about half of L(< rcool)- 



Edge (1987) and White et al. (1997), but to our knowledge, 
it is the first systematic investigation on the fraction of a 
cluster luminosity contributed by its core, undertaken with 
ROSAT. 

We note that the large contribution to the X-ray bolo- 
metric luminosity emerging from within the inner -^200 kpc 
of a cooling flow cluster must have an influence on its de- 
tectability in imaging surveys (Pesce et al. 1990) and on the 
scatter in the Lx — Tx correlation (Fabian et al. 1994) . 



4.2 New Cooling Flows in the Sample 

We present in Table ^ the mass deposition rate, cooling ra- 
dius, and central cooling time for the clusters which did not 
feature this information in the Edge et al. (1992) paper. 
Three interesting cases are A2065, A1650 and A2204. In the 
first two cases the resolution of the data in the Edge et al. 
paper hid a modest and a large cooling flow. The third case 
(A2204) is a massive cooling flow. A summary of the depro- 
jection results for this cluster is presented in Fig. H (A2204 
was observed with the ROSAT PSPC from 1992 September 
04 to 1992 September 05 and with the ROSAT HRI from 
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Figure 3. The luminosity from the cooUng region, L{< rcool)) 
is plotted against the X-ray bolometric luminosity from David et 
al. (1992). We note that, despite the dispersion, many clusters 
have L{< r^^^l) contributing to 50 percent or more to the total 
X-ray luminosity. The dotted lines represent 20, 50 and 100 per- 
cent of the cluster X-ray bolometric luminosity. Arrows represent 
non-cooling-fiow clusters which were assigned L{< r;.ool)=0.01 
for display purposes only. 



Mass Deposition Rate (Mq yr ) 

Figure 4. /loo is defined as the luminosity from a region of 
100 kpc radius divided by the cluster X-ray bolometric luminos- 
ity given by David et al. (1992). We note that non-cooling-flows 
(displayed here with M = IMq yr~^) have /loo < 10 percent, 
whereas cooling flow clusters have values of /loo ranging from 5 
to 50 percent. A trend of increasing /loo with mass deposition 
rate is indicated by the data. 



Table 4. New cooling flows in the B55 Sample. The values 
quoted for M, r^ool^ a-nd tcool a^re median, 10 and 90 percentile 
estimates from the deprojection algorithm. The units used for M, 
''cooli and tcool are WIq yr"'^, kpc, and Gyr respectively. 

Cluster z M r^ooi tcooi 



Klem44 


0.0283 


871^5 


133t« 


r, O-I-0.6 

^■•^-0.3 


A1650 


0.0845 




165j:if 


2 4+1-2 


A1651 


0.0846 




127^31 


fi 5+°-^ 


A2204 


0.1523 






o 1 -t-0.1 

■'■^-o.i 


A2065 


0.0722 




561^3 


4 4+2-2 


A3558 


0.0478 




681^^ 


10 2+°-^ 



1995 January 01 to 1995 January 11, with total exposure 
times of 5.36 ks and 15.5 ks respectively.) 

At the same time that new cooling flows were discov- 
ered three clusters considered to have small flows are now 
classified as non cooling flows: A2319, A576, and A754. 
A2319 does not have the single peak typical of a cooling 
flow when observed with the HRI. The central cooling time 
is marginally consistent with the assumed age, and it poss- 
eses a radio halo (Hanisch 1982). A754 has a very disturbed 
X-ray morphology and a central cooling time in excess of 13 
Gyr. It is known to be in a merging state (Henry & Briel 
1995, Zabludoff & Zaritsky 1995, and Henriksen & Marke- 
vitch 1996). A576 has a very flat surface brightness in X-rays 
when observed with the HRI. Although the cooling time is 



low the mass deposition profile is perturbed and the value of 
the mass deposition rate is consistent with zero. From our 
analysis a clear cooling-flow/morphology relation appears 
in the B55 sample. This is consistent with the results of 
Buote and Tsai (1996) who found a quantitative correlation 
between the mass flow rate and morphology of 37 ROSAT 
PSPC clusters of the B55 sample. A similar analysis for the 
whole B55 sample will appear elsewhere (Peres, Buote & 
Fabian, in preparation.) 

All clusters which we classify as non-cooling-flows dis- 
play signs of merger activity and/or a disturbed morphology 
in X-rays. Most notably, there are no examples, in the B55 
sample, of a cooling flow cluster without a well defined cen- 
tral galaxy at the bottom of the cluster potential well 
The lack of ROSAT data for A2147 made us keep its clas- 
sification as a non-cooling-flow cluster, since it had a large 
central cooling time and a mass deposition rate consistent 
with zero in the analyses of Edge et al. (1992) and White et 
al. (1997). 

4.3 Fraction of Cooling Flows 

The work of Edge et al. (1992) addressed this issue within 
the B55 sample, concluding that the fraction of cooling flow 
clusters is 70 percent. However they warned that it could be 
"90 percent, once the effect of different bin sizes is taken into 

** The existence of non-cooling-flow clusters with a cD galaxy at 
the bottom of the potential well, e.g. A399 and A401, is very rare 
leading to the interpretation of it as evidence of a previous phe- 
nomenum, which disrupted the flow (McGlynn & Fabian 1984; see 
Fabian, Peres & White 1997 for the specific case of A399/A401). 
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Figure 5. Summajry of the deprojection results for the PSPC observation of the cluster A2204. Points in Figs (a), (b), (c), (e) and (g) 
represent the mean value and Icr errors (in each radial bin) from 100 Monte-Caxlo simulations. Points in Figs (d), (f) and (h) represent 
the median and 10 and 90 percentile estimations from 100 Monte-Carlo simulations. The dotted line in panel (d) marks the average 
temperature from David et al. (1992), and Rcen denotes the radius at the bin centre. Note that the potential is adjusted such as to 
produce a flat temperature profile. 
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Figure 6. Central cooling time plotted against binsize for the 
B55 sample. The dotted line indicates the assumed age of the 
Universe (13 Gyr). From this plot it is evident that 70-90 per- 
cent of the clusters in this sample are cooling flows. The trend 
of reduced central cooling time with increased resolution is also 
easily observed for the sample as a whole. The stars and squares 
represent PSPC and HRI observations respectively. 



accounf . In the study by Edge et al. (1992) only 36 clusters 
(65 percent) had observations with a resolution better than 
100 kpc. Correction for that was made on the basis of a 
tendency for the cooling time to decrease with bin size as 
tc ~ R}'^ , which the authors inferred from their sample. 

Here we confirm this tendency (cf. Fig. but do not 
resort to it to establish the fraction of cooling flows in the 
B55 sample. With the improved spatial resolution of ROSAT 
we can resolve the central cooling regions of all clusters in 
the B55 sample as can be seen from the short cooling times 
displayed in the histogram in Fig. ^. 

We define here a cluster as a non-cooling-fiow when the 
upper limit (90th percentile) to the central cooling time, as 
given by the deprojection procedure, is larger than the as- 
sumed age of the cluster (cf. Sec. ^ and Figs. ^ ^ above). By 
doing so we conclude that the fraction of cooling fiows in the 
B55 sample is 70 percent (consistent with the 70-90 percent 
estimates of Edge et al. 1992). Importantly, this is a con- 
servative value since we are not considering as cooling fiows 
clusters for which the lower limit to the central cooling time 
is consistent with their assumed age (i.e. 13 Gyr). If we did 
include the clusters for which the central cooling time is con- 
sistent with 13 Gyr then the fraction of fiows would be ~90 
percent and we would have to include as cooling-fiows A119, 
A401, A399, A3158, A3266, A754, A3532, A1736, Triangu- 
lum Australis, A2256, A2319, and A3667. We prefer the first 
possibility here because it does not include as cooling fiow 
those clusters which are clear mergers, but whose central 
cooling time is just above the assumed age of the Universe; 
the ages of these clusters are certainly less than the latter. 

In hierarchical scenarios for the formation of structures 
in the Universe clusters are assembled from smaller subunits 
through mergers. This complicates the estimation of the age 
of a cluster (and of its cooling flow age) , which is one of the 
parameters used to determine the mass deposition rate. The 



fact that a cluster has suffered a merger does not necessarily 
imply that its cooling flow has an age equal to the lookback 
time of the merger, i.e. cooling flows may be older than their 
recently formed host cluster; as yet we do not know in de- 
tail how a merger affects an existing cooling flow |^. In fact 
if only approximately equal-mass, head-on collisions disrupt 
existing cooling fiows then the rarity of such events is recon- 
ciled with the high fraction of flows in our sample and can 
be used in future larger datasets to constrain cosmological 
scenarios. 

A detailed investigation of the sensitivity of the mass 
deposition rate and fraction of cooling flows to individual 
cluster age would be ideal, but the complex history of for- 
mation of a cluster justifies the use of an average cluster 
age in our analysis. Assuming an average age of 13 Gyr for 
all clusters we obtained a fraction of cooling flows of 70-90 
percent in the B55 sample. Since clusters are younger than 
13 Gyr in all plausible formation scenarios M and the frac- 
tion of cooling flows derived here are upper limits to their 
true value. Assuming an average cluster age of 6 Gyr yields 
a fraction of cooling flows equal to 65 percent, whereas an 
age of 2 Gyr would still classify 45 percent of our clusters as 
cooling fiows. The values of M do not change by more than 
a factor of three when we assume an average cluster age of 
5 Gyr (With the exception of A644, A1650, A1651, A2142 
and A4059, they either change by less than a factor of two 
or do not change appreciably.) 

The high fraction of cooling fiows found here is ob- 
jectively defined and secure. The sample is not biased by 
the presence of large cooling fiows because it is fiux-limited 
and constructed from observations with broad-beam instru- 
ments. All the clusters considered as cooling fiows have very 
short central cooling times, and the fraction of flows in the 
sample remains basically unaltered if any other (reasonable) 
cosmology is selected {tc/Hg^ scales as Hq'^; cf. Fabian et 
al 1984). 

We have also computed the cooling time of the hot 
phase at a fixed radius (250 kpc) for as many clusters in 
our sample as possible. The value quoted in Table ^ t250, 
is the cooling time of the bin encompassing the 250 kpc ra- 
dius and the quoted limits come from the Monte-Carlo im- 
plementation of the deprojection analysis. The results are 
summarised in Fig. ^, from which is clear that all but two 
clusters have cooling times 13 < icooi(250 kpc) < 51 Gyr. 
Fig. ^ was produced under the approximation that the lim- 
its to t25o were la extrema of a gaussian distribution (the 
limits were computed from the uncertainties provided by the 
deprojection analysis). 

This exemplifies once more the need for high spatial res- 
olution in the determination of a cooling fiow and provides 
a first estimate of the range of cooling times at around 250 
kpc required from future numerical simulations of clusters 
of galaxies. The short cooling times at 250 kpc found in our 
sample are at odds with the values in the simulations by 
Cen et al. (1995) and support the multiphase picture of the 
ICM where there can be gas cooling at and beyond 250 kpc. 
Future higher resolution data may demonstrate that some 
of the the B55 clusters which have central cooling times 



tt A85 is the example of a cluster with clear substructure which 
harbours a strong cooling flow. 
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Figure 7. Cooling time at 250 kpc plotted against cumulative 
number of clusters. The plot shows 48 clusters for which the 
deprojection analysis was undertaken to a radius greater than 
250 kpc. Each cluster was represented by a normalized gaussian 
with dispersion given approximatelly by the deprojection limits. 
This plot shows that cooling times at 250 kpc from the centre of 
cooling-flow clusters are higher than the age of the Universe (13 
Gyr) by not more than a factor of ~4. 



larger than their assumed agee (13 Gyr here) are nonethe- 
less cooling-flows, as indicated from the trend in Fig. ^ (al- 
though the trend is not guaranteed for every cluster individ- 
ually, meaning that there are clusters for which we do not 
expected to detect a cooling flow even with data of improved 
resolution) . 



4.4 Mass Deposition in Cooling Flows 

Cooling flows are known to deposit large amounts of cool 
material throughout a cluster core (cf. Sec. i.i. and Fig. |^. 
From table |5| is clear that the accreted mass during the age 
of the flow. Mace ~ M X tagc, varies from 6 x 10^° M© to 
6 X 10^^ Mq , if we assume an average cluster age of 6 Gyr and 
meiss deposition rates of 10 and 1000 M© yr~^ respectively. 
We note from the discussion in the previous section that 
these estimates are not altered by more than a factor of 
^ 6 when the cluster age is halved since M does not change 
by more than a factor of ~ 3. The mass accreted in large 
cooling flows (> 2OOM0yr~^ ) should not change by more 
than ~ 2. Therefore it would be desirable to understand how 
this material is distributed over the cooling region, and how 
the cooling cores are distributed in the Universe. 

We have studied the spatial variation of this material 
over the cooling regions of the B55 clusters with the HRI. To 
avoid contamination by small scale substructure we use the 
integrated mass deposition (IMD) proflle which comes from 
the deprojection analysis and therefore has information on 
the azimuthally averaged mass deposition only. The shape 
of the IMD profile can be used to test the multiphase nature 
of the IGM and to trace the evolutionary history of the flow. 

As for the flrst case, we know from previous work 
(Nulsen 1987, Thomas et al. 1987) that simple models of 
a multiphase flow would lead to an IMD proflle of the form 
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Figure 8. Histograms of mass deposition rates (top) and central 
cooling times (bottom) for the ROSAT observations of the B55 
clusters. Note that changing the value of Hq does not affect the 
fraction of cooling flows in the sample, because the ratio tc/H^^ 
scales as i/Jj"'' (cf. Fabian ot al 1984). 



A/(< r) ~ r" , with a ~ 1. We have used HRI and some 
PSPC observations to compute the slope of the IMD pro- 
file. For this purpose we have obtained the profiles from the 
deprojection algorithm with data points given by the mean 
and la deviations. The profiles were then fit with a power 
law model. Only clusters with 4 or more bins inside the 
cooling radius were considered in the fitting procedure. Av- 
eraging the results for the whole set of observation analyzed 
here yields (a) = 0.95 and a histogram of the values ob- 
tained can be seen in Fig. ^. The data is in good agreement 
with the simple multiphase models for the ICM. (Arnaud 
1987 reached similar conclusions from a different sample.) 

The possibility of using the IMD profiles as tracers of 
the evolutionary history of the flows was raised by Edge et 
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Figure 9. Histogram with the values of the slope of the integrated 
mass deposition profile for the B55 clusters. Note that the values 
scatter around the value a ^ 1 predicted by multiphase models . 



al. (1992). They searched for the existence of 'breaks' or 
'plateaus' in the IMD profiles, as indicative of previous phe- 
nomena which perturbed the fiow. The detection of such 
features in the IMD profiles depends on how well we can 
sample the cooling region, for which HRI data is ideal. To 
search for a 'break' in the IMD profiles we fit them with a 
broken power law in the same manner described above. We 
have noticed that many clusters show some kind of break in 
their IMD profiles but only few are as apparent as in Fig. |l^ 
(see also Irwin & Sarazin, 1996); the clusters showing clear 
breaks in our analysis are 2A0335+096, A426, and A496. 
In the case of 2A 0335+096 the coofing time at around 70 
kpc, where the break occurs, is ^ 4-5 Gyr; in the sugges- 
tion by Edge et al. (1992) this is the approximate age of 
the flow (time since last disruption). Some of the breaks, 
however happen in the outer bins of the IMD profiles, which 
can be attributed to fluctuations in the temperature proflle, 
since M goes to first order as Lx/T. Allen and Fabian (1997) 
have proposed a method to estimate the age of cooling flows 
based on the spectral capabilities of the PSPC which needs, 
however, high signal-to-noise data. Future high spatial res- 
olution missions like AXAF will pin-down this issue. ^ 

We have also used our recent data on the mass depo- 
sition by cooling flows to recalculate the space density of 
the flows (Edge et al. 1992). We have selected the clusters 
with galactic latitude larger than 20 degrees and applied 
the usual procedure to compute the density function for a 
flux-limited sample. The result is presented in Fig. nil; the 
slope of a power law fltted between M = 10 M© yr~ and 
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-1.8. Selection of clusters with 



Numerical experiments are currently being undertaken with 
the aim of understanding how cluster collisions can affect the 
central cooling flow (e.g. Burns et al. 1994, Roettiger et al. 1997). 
Once the cause of disruptio n o f a cooling flow is known, the inci- 
dence of cooling flows (sec. [4.3| ) may be used to constrain cosmo- 
logical scenarios. 
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Figure 10. Integrated mass deposition (IMD) profile from the 
deprojection of an HRI observation of 2A 0335-1-096. Note that 
the slope changes clearly inside the cooling radius for the cluster 
(fcool ^ 200 kpc). Error bars are 10th and 90th percentiles from 
our deprojection analysis. 
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Figure 11. Cumulative space densitive of cooling flows from the 
B55 sample. The density is corrected for the selection function 
|6| > 20 degrees. The space density of non-cooling-flows, which 
were assigned here a mass deposition rate equal to one, is repre- 
sented by a star. 



\h\ > 10, instead of 
significantly. 



6| > 20 does not change the results 



4.5 Gas fractions in cluster cores 

Gas masses within 500 kpc or 250 kpc for the B55 sample are 
presented in Table |^ They range from 2.5 to 64x10^^ Mq, 
constituting from 7 to 23 percent of the total gravitational 
mass within this region. 

Systematic errors due to our choice of potential and 
methodology are larger than the limits quoted from the de- 
projection analysis. However, since the baryon fraction in 
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Gas fraction (percentage) 

Figure 12. The histograms above show the baryon fractions for 
non-coohng-fiow clusters (top) and for coohng-fiow clusters (bot- 
tom). The dotted line in both panels represents the cases for which 
baryon fractions were calculated within 250 kpc only, whereas the 
solid line includes clusters with gas fractions computed within 500 
kpc. 



clusters (/b) tend to increase with radius, we note that for 
our sample, /b > 6 percent. This conclusion is robust to the 
systematics mentioned above, agreeing with previously pub- 
lished resuhs (White et al. 1993, and White & Fabian 1995). 
It is important to stress that the gas fractions presented here 
do not include the contribution of baryonic material in clus- 
ter galaxies, representing lower limits to the baryon fractions 
in clusters. 

We plot in Fig. ^ the gas fraction for cooling-flows and 
non-cooling-flows within 500 kpc from our analysis (250 kpc 
when information for 500 kpc is not available). This plot 
indicates that the non-cooling-flows have a similar amount of 
gas, for a given mass, to cooling-flow clusters; the difference 
between them residing therefore in the spatial distribution 
of the gas in the core. 



previous analysis it becomes clear that: (i) cooling flows are 
common and long-lived (fraction of at least 70 percent in a 
flux- limited sample), and (ii) they deposit a large amount of 
cool material throughout a cluster core (about 40 percent of 
the clusters in the B55 sample have M > 100 Mq yr~^ ). 

As a corollary of these points we state that the cooling 
flow plays an important role on practically all phenomena 
happening in a cluster core. Historically, two phenomena 
have been traditionally linked to the cooling flow: optical 
emission line nebulosity and radio activity from the BCG. 
A review of the role of cooling flows in these phenomena can 
be found in Fabian (1994). 

Here we use the data from the work of Heckman et 
al. (1989) to investigate the role of the flow in the Ha lu- 
minosities of the BCG. We have changed the values of the 
published luminosities to match our assumed cosmological 
parameters. The results of a cross-correlation between the 
B55 and Heckman's sample (24 objects in common) is dis- 
played in Fig. |l^. We note that the weak trend suggested 
by Heckman et al. is maintained, i.e. Ha luminous galaxies 
lie in the centre of large cooling flows although this spe- 
cial cluster environment does not guarantee the existence of 
emission line nebulosity in its BCG. Finally we note from 
Fig. |l4|that clusters which display emission line nebulosities 
have the shortest central cooling times in the sample: only 
two greater than 2 Gyr. 

Various generic scenarios have been put forward to un- 
derstand the trend in Fig. ^ Crawford & Fabian (1992) 
suggest that the emission line nebulosity is powered by 
mixing-layers in the multiphase, turbulent ICM (Begelman 
& Fabian, 1990). Since mixing layers occur where there is 
a population of cold clouds and where the hot ICM is most 
turbulent, we would expect emission-line nebulae to occur 
in the centres of large cooling flows. On the other hand this 
scenario is consistent with the observation that some large 
cooling flows do not host emission-line nebulae (e.g., A2029) 
since the hot ICM may be in a quiescent state and/or most of 
the cold clouds may have formed stars already. It is interest- 
ing to remember that A2029 shows little excess absorption 
(White et al. 1991, Allen & Fabian 1997). Allen (1995) has 
presented new evidence for star formation at the centres of 
cooling flows, proposing that it can power the emission line 
nebulae; early suggestions of this mechanism are found in 
Johnstone et al. (1987), Hu et al. (1985), and Heckman et 
al. (1989) among others. The role of magnetic reconection on 
the emission line nebulae have been re-investigated recently 
by Jafelice & Friaca (1996), and a recent detailed case study 
has been presented by Voit & Donahue (1997). 

Optical images retrieved from the Space Telescope Digi- 
tized Sky Survey (DSS) were used to locate the optical coun- 
terparts to the peak in the X-ray emission from our cluster 
images. The optical position listed in Table ^ is the opti- 
cal peak of the brightest galaxy nearest to the position of 
the X-ray peak. The listed offsets, AS, are the modulus of 
the difference between the X-ray and optical positions, but 



4.6 Comparison to Optical and Radio Data 

Cooling flows are the only simple route to explain the wealth 
of X-ray data on the centres of clusters of galaxies. From our 



38 About 50 percent of the cooling flow clusters in the sample 
show line emission, whereas about 75 percent of the clusters with 
emission lines display cooling flows with M > 100 Mq yr~^ . This 
trend can be used as an efficient method to search for large cooling 
flows at high z. 
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no effort was made to correct the astrometry of the X-ray 
images. We note that most of the clusters with large off- 
sets, AO > 25 arcsec, do not have cooling flows: A119, A401, 
A3158, A754, A1367, A1736, A2256, and A3667. Most of the 
clusters, however have optical positions consistent with the 
uncertainty in the position of their X-ray peaks: 65 percent 
have A9 < 10 arcsec, and 85 percent have A9 < 25 arcsec. 
The large values of A9 for A1367 and A1736 stem from the 
flat X-ray surface brightness, without a clear centre. The 
extreme offset values of A2256 and A754 are caused because 
the peak of the X-ray brightness does not correspond to the 
BCG in these merger systems, but to shocked regions of the 
ICM. The overall excellent agreement between the peak of 
the X-ray emission and the position of the BCG indicates 
that for most of the clusters in the B55 sample the X-ray 
gas is in a dynamically quiescent state, with the centre of 
the potential being occupied by a bright galaxy. 

The appearance in the literature of a new, larger, pho- 
tometrically homogeneous sample of BCGs (Lauer & Post- 
man 1994) has prompted us to use our new cooling flow 
results to search for a correlation between the strength of 
the flow and the optical light from the BCGs. This question 
has been investigated in the past in relation to the existence 
of a putative accretion population formed from the flow by 
Schombert (1987), Sarazin (1983), Thuan & Puschell (1989), 
and Edge (1991). 

The sample by Lauer & Postman is volume-limited to 
cz <15,000 kms"\ and it has 19 BCGs (out of 119) in 
common with our sample. We plot M against the absolute 
magnitude in Fig. from which we notice that no evident 
correlation exists. The lack of correlation can be understood 
if star formation is skewed towards low-mass objects in the 
environment of cooling- flows, but this is a much harder ques- 
tion to answer observationally (cf. Fabian 1994). 

We have also cross-correlated the strength of the flow 
with the radio power from BCGs. We have used radio data 
on 17 BCGs by Ball et al. (1993) and archival data obtained 
from NED, and the NVSS survey. From the archival informa- 
tion we could obtain the flux density of the central source 
and then its luminosity at 6-cm and/or 20-cm. No simple 
correlation is evident from Fig. |l^ (see also Fabian 1994), 
but as we see, large cooling flows tend to have high radio 
luminosities (right upper corner of the diagram) . This trend 
can be understood as a consequence of the high-pressure en- 
vironment of cooling flows, which provide the working sur- 
face needed to produce the radio emission. When we split the 
data between clusters showing and clusters lacking emission 
lines, we note that the trend exists for the former only. In- 
spection of Table ^ shows that 19 out of the 22 clusters (i.e. 
90 percent) which display emission lines have radio emission 
from the center, although only 19 out of 27 clusters (i.e. 70 
percent) displaying central radio emission have emission line 
nebulosity. 

Finally we note a similarity between Figs. ^ and |l6[ 
Both have a distinctive upper envelope with a roughly uni- 
form distribution of data under it. In neither figure do we 
see large cooling flows with a central cluster galaxy of small 
absolute magnitude or small radio power. 
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Figure 13. Correlation of Ha luminosities from the centres of 
clusters listed in Hcckman et al. (1989) and tlie mass deposition 
rates inferred from our deprojection analysis. Note that for clarity 
we did not include the Virgo cluster in this sample. Points with 
left-pointing horizontal arrows represent upper limits for the ob- 
served Ha luminosity. 
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Figure 14. In the diagram above we denote clusters with emis- 
sion lines by squares and clusters lacking emission lines by filled 
stars. Open squares are used to represent the clusters for which we 
have no information on the presence of emission lines. Emission 
line clusters are situated predominantly below the line of central 
cooling time equal to 2 Gyr. Note there are only two emission line 
cluster with central cooling times larger than 2 Gyr. 



5 CONCLUSIONS 

Our analysis of the X-ray properties of cluster cores with 
ROSAT showed that a cooling flow is the natural state of 
the cores of nearby clusters of galaxies. They constitute 70- 
90 percent of the clusters in our flux-limited sample and 
deposit more than 100 M0yr~^in about 40 percent of the 
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Figure 15. The lack of correlation displayed in the figure above 
indicates that the putative accretion population from the flow is 
not composed of stars with a normal IMF. This however does 

not eliminate the possibility of ongoing low-mass star formation. 
Despite the lack of correlation, there is the indication of a cut-off 
in the upper left part of the plot. Non-cooling-flows were assigned 
M = 1 for display reasons only. 
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Figure 16. The plot above shows a correlation between the inte- 
grated power at 6cm emitted by a cluster (Pgcm) and its cooling 
flow strength. Points to the right of the horizontal arrows repre- 
sent upper limits for the observed radio power. Stars represent 
clusters with observed optical line emission; squares are used to 
denote absence of emission lines or lack of information thereof. 
Non-cooling-flows were assigned M = 1. Data for which only up- 
per limits exist were not used in the plot. 



clusters. Large cooling flows can contribute more than 70 
percent of the cluster bolometric luminosity. 

Prom our catalogue of ICM properties of the B55 sam- 
ple we note that cooling times at the central bin are small 
and cooling times at 250 kpc from the cluster centre do not 
exceed ~4 times the age of the Universe (13 Gyr). We show 
that mass deposition in cooling flows have an approximatelly 
linear dependence on radius and that clear breaks in some 
of the integrated mass deposition (IMD) profiles exist. 



The cross-correlation of optical and radio properties 
with the strength of the flow yielded weak trends only. We 
noted that overall large cooling flows tend to have luminous 
BCGs and powerful radio sources at their centres, whereas 
no such rule exists for non-cooling-flows. 
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Table 5. Summary of some deprojection results, (a) Name of the cluster. The letter inside brackets indicates wether the observation 
was made with the PSPC(P), or with the HRI(H); (b) Integrated mass deposition rate in units of Mq yr~^ ; (c) Cooling radius in units 
of kpc; (d) Central cooling time in units of Gyr; (e) X-ray bolometric luminosity within the cooling radius, L(< r^ooi), in units of 
10^^ ergs~^ ; (f) Ratio of L{< rcool) to Lbol- The values of Lbol s^re from David et al. (1992); (g) Cooling time of the bin encompassing 
the 250 kpc radius, in units of Gyr; (h) Gas mass out to 0.5 Mpc or 0.25{t) Mpc, in units of 10^"^ Mq! (i) Gravitational mass out to 0.5 
Mpc or 0.25(t) Mpc, in units of 10^^ M©; (j) Approximate ratio (percentage) of gas-to-gravitational mass inside 0.5 or 0.25(^) Mpc. 
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A2029(P) 




1861^^ 


2 9+°-^ 




42+0 03 


2i-5^i;? 


44.5±0.5 


251 


18±0.2 


A2052(P) 


125tf 


147^3^ 


2 5+''-l 

^■°-o.i 


-^■''-0.02 




26.0t-« 


15.9±0.4 


107 


14±0.4 


A2052(H) 


102tif 




1 q+0.3 

^■'-'-0.2 


1 1+0.2 
^•'-0.1 


0.46t«j 


32.01:^571 


5.1±0.3 


49 


10±0.6 
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Table ||- cont. 



Cluster 


M 




^cool 


»'cool) 


iCO-cool) 


i250 




-^^grav 


Mgas 


^cool 






(a) 


(b) 


(c) 


(d) 


(c) 


(f) 


(g) 


(h) 


(i) 


(j) 


MKW3(P) 


. 4_1 A 


1 yi -t-y 

-^'-^-62 


„ ^-Lfl 1 

Q n ' 
■^•^-o.i 


O o + O.l 

^■"^-0.6 


„ „„-i-n (19 

0.52+0-02 


22.7+1-g 


7.1±0.lt 


58+ 


12±0.2t 


MKW3{H) 


LI r;n 


137+7 


_ _i_n Q 
^•^-0.3 


.4-0 d- 

^■^-0.5 


0.32+O;^j5' 










A2065(H) 


13^6 




. A 4-2 2 
^•^-1.3 


„ „4-n 9 

0.5l°-J 


^ ^ . 4~n 01 
"-"•^^-coi 


*'-"-16.0 


6.8±0.5t 


90+ 


8±0.6+ 


A2063(P) 




^■^ -1-1 

95t35 


^ ^-1-0 4 

5-oi;j:3 


„ ^4-0 1 
'^■'^-0.3 


„ „„-i-n 0*^ 
0.20+0 J 


30-8-2:9 


15.8±0.3 


116 


14±0.3 


A2142(P) 






c- fj + 0.4 
^•^-0.3 


^'^■'-'-4.6 


n 24+0-03 


1 o + l^-8 
^'•^-0.8 


52.0±0.8 


369 


14±0.2 


A2142(H) 


286+^4 


152+23 


„ n + 1 9 

3-8_0'8 


, _ „4-9 1 




^-l-fi s 
26.61^2 


47.2±0.3 


349 


14±0.1 


A2199(H) 




_ „ „ _!_ 1 7 


1 9t^'-^'-== 
-^•^-0.1 


„ „4_n 9 
■^■^-0.5 


^•^'-0.08 










A2199(P) 


_ „ . -1-1 s 

154Ig^ 


, . 7 

143Ig 


^ „-i-n n 
-^•^-0.1 


-|-0_3 

^- ' -o'.i 


^•^^-0.02 


^4-1 3 

27-9_i:5 


20.8±0.2 


146 


14±0.1 


A2204(P) 


n-f, 4-127 

852+if 


199^44 


o -1 4-0 1 
-^•^-0.1 


, „ ^ 4-4 

40- Ills 


r\ --T----)-0 09 

0-75Io.o9 


^''•^-1.6 


51.6±1.2 


339 


15±0.4 


A2204(H) 


843t2« 


181«« 




40.3^4 2 




23-8j:I:5 


21.0±0.6t 


186t 


11±0.3+ 


TRI AUST(P) 


QQ-t-O ( 

■^•^-33 


-Lfi7 


-L7 1 

io.8i^:i 


-LI 7 


-Ln rifi 

0.03IoO6 


4-1 Q 

24.411;^ 


38.8±0.7 


197 


20±0.4 


A2244(P) 




148^20 


7.3i»:« 


"■"-4.0 


^ . ^4-0 0"^ 
'-'•^'^-0.26 


27.91:3-2 


35.7±1.2 


284 


13±0.4 


A2256(P) 




0+ 




in 9 

0+ 
^-0 


„4_n ni 


35.5+|;g 


24.4±0.4 


154 


16±0.3 


A2256(H) 


(1+16 


ii;+69 
^"-16 


10.8tg°i^ 


O.oto^" 


n+o.o 

^-0.0 


28.015-2 


26.5±0.3 


159 


17±0.3 


OPHI(P) 


1 27+ 


1 29+ 


4 

■^■'^-0.3 


4-0 Q 

Q Q-|-u.y 
3.»_2 g 




4-1 Q 

20.lll;«g 


37.6±0.7 


242 


16±0.3 


OPHI(H) 


-Lin 
341 J° 


-L7 

72Ii 


1 

1 n+ 


-Ln 1 
-^•'-'-0.01 


0.04 










A2255(P) 


U_o 




^-4-91 n 
36-9tg'5 


^-0 


n+0.01 


0+7.0 

54-2-6:9 


21.5±0.5 


171 


13±0.3 


A2319(H) 


-1+103 


^ on 
34t™ 


„ . -LI 

^•^-3.7 


^ -. 4-9 fi 

0.1+p" 




si-el^'e' 


9.9±0.5t 


103+ 


10±0.5+ 


A2319(P) 




53tf3 




0.51J-5 


O.Ol^oo* 


26.611° 


37.0±1.2 


240 


15±0.5 


CYG-A(H) 


2151^6 




-Ln 1 
-^•'-0.1 




-LO 09 

^•^^-0.02 


4-Q A 

28.5l^;^ 


9.0±0.3t 


109+ 


8±0.3+ 


CYG-A(P) 


„ . . 4-26 

244+^^ 


135+33 


r. ^4-0 1 

^••^-0.1 


„ „_Ln 8 

9.31^-" 


^ _„4-0 05 

0.56+0-05 


— L 3 fl 

28. 5+21 


28.2±0.8 


201 


14±0.4 


A3667(P) 






^•^•^-3.4 


4-0 9 
^-0 


'-'-0 


4-9 Q 

28-7l2:3 


24.8±0.4 


153 


16±0.3 


A2597(H) 


276_3 


139^24 


-^•^-0.1 


8-2lo.9 


„ „ ^ -1-0 
'-'•°^-0.07 










A2597(P) 


271+« 


152^«^ 


^•■^-0.1 


8-2_i.9 


61+0-1 


32.7l2;3 


22.9±0.6 


217 


11±0.3 


KLEM44(P) 




133t^? 


r, O + 0.6 

^•■^-0.3 


1 1+0.3 
-^■-^-0.2 


49+0-13 


00 0+121.0 

»3.S_42 7 


4.3±0.2t 


50+ 


9±0.4+ 


A4059(H) 






1 0+1-14 


1 5+°--* 


r, 09+0.09 
"•'-'^-0.02 










A4059(P) 


130+2^ 


153^^3 


3 4+0.6 


1 6+°-^ 


34+0-04 
"•'-'^-0.04 


24.8ll« 


19.6±0.6 


98 


20±0.6 



Notes. The errors quoted for rcool E're smaller than the resolution used. This must be understood remembering that they are 
uncertainties computed from the deprojection of Monte-Carlo perturbations of the X-ray surface brightness. The limits in L{< r^aol) 
correspond to the bolometric luminosities within the extrema of rcool- We use the integrated X-ray luminosity at the last deprojected 
bin as the X-ray bolometric luminosity for A3532, since the latter is not available in David et al. (1992). The luminosity for M87 in the 
Virgo cluster is taken from White et al. (1997). 
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